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Abstract. Wc discuss the properties of the net-quark and isovector fluctuations along the chiral phase 
transition line in the plane spanned by temperature and baryon chemical potential. Our results are obtained 
in terms of the Nambu— Jona-Lasinio (NJL) model within the mean-field approximation. The model is 
formulated at the finite temperature and for non- vanishing net-quark and the isospin chemical potentials. 
The fcrmion interactions are controlled by the strength of the scalar and vector couplings in the iso-scalar 
and iso-vector channels of constituent quarks. We explore properties and differences in the behavior of the 
net-quark number and isovector susceptibilities for different values of thermal parameters near the phase 
transition. We argue that any non-monotonic behavior of the net— quark number susceptibility along the 
phase transition boundary is an excellent probe of the existence and the position of the second order endpoint 
in the QCD phase diagram. 

One of the essential predictions of QCD is the existence of the boundary line in the temperature and 
net-quark chemical potential, (T, /^ g )-plane that separates the confined, chirally broken hadronic phase from 
the deconfined quark-gluon plasma phase. The existence of such a boundary line for fX q /T < 1 has been 
recently established by the first principle calculations in the Lattice Gauge Theory (LGT) formulated at 
finite baryon density P El E] . 

Arguments based on effective models El El EH QI3 ^2 U2\ Q3| indicate that at large \i q the transition 
along the boundary line is the first order. For small fi q and for two masslcss flavor QCD the chiral transition 
was argued ^3] to be second order with the critical exponents of the 0(4) spin model. For the finite quark 
mass, due to the explicit chiral symmetry breaking, this transition is likely to be replaced by the rapid 
crossover. Such a different nature of the phase transition at low and high fj. q suggests that the QCD phase 
diagram should exhibit a critical endpoint at which the line of the first order phase transition matches that 
of the second order or analytical crossover The critical properties of this second order chiral endpoint 
are expected to be determined by the three-dimensional fsing model universality class llfij . 

The existence of a critical endpoint in QCD has been recently studied in the lattice calculations at the 
non-vanishing chemical potential by cither considering the location of Lee- Yang zeros in (2+l)-flavor QCD 
I17j or by analyzing the convergence radius of the Taylor series in fi q /T expansion of the free energy 
1181 HT?j . Recent results pQ based on the first approach suggest that a critical endpoint indeed exists in 
QCD phase diagram and might occur at T ~ 164 and (i q ~ 120 McV. In the 2-flavor QCD and with a 
relatively large quark mass used in the actual lattice calculations ^Hj no direct evidence for the existence of 
the critical endpoint has been found for \x q < T where the Taylor expansion method is applicable. 

The critical behavior and the position of the chiral endpoint can be possibly identified by observables 
that are sensitive to singular part of the free energy |20| . One such observable is the quark susceptibility Xij 
defined as the second order derivative of the thermodynamical-potential 0(T, fi, V) with respect to quark- 
flavor chemical potentials, Xff = —d 2 Q,/ d^ifd^f , where for two light (u,d)-quarks, /i = {[i u ,Hd)- To 
identify the chiral critical endpoint in the QCD phase diagram, the properties of the net-quark number 
susceptibility \ q are of particular interest |21l 1221 12"3| . In two-flavor QCD, \q is expressed as the sum of u 
and d quarks susceptibilities: \ q = 2{Xuu + Xud)- The universality argument predicts that independent of 
the values of the quark masses the Xq should diverge at the chiral endpoint. 

The quark susceptibilities have been recently obtained ^1 11DI within the lattice QCD for two 
light quark flavors using p4-improved staggered fermion action with the quark mass m q /T = 0.4 0EI1- The 
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Figure 1. The left-hand figure shows the LGT results obtained in 2-flavor QCD 1181 1251 for net-quark 
number susceptibility \q as a function of the temperature normalized to To being the transition temperature 
at y,q = . The results are shown for different values of /J,q/T. The dashed-dotted lines indicate the 
corresponding hadron resonance gas model results 1251 . The right-hand figure: as in the left-hand figure but 
for the isovector susceptibility. 



Monte Carlo simulations have been done at the finite quark chemical potential via Taylor series expansion. 
The susceptibilities were calculated up to the 0([i q ) order in the quark chemical potential ^0 120] ■ 
The lattice results ^0 for the net-quark number \q an d isovector susceptibility xi as a function of 
the temperature for different values of [i q are shown in Fig. 1. There is a strong suppression of the Xq 
fluctuations in confined phase and the cusp-like structure in the near vicinity of the transition temperature 
To- A rather strong increase of Xq a t the transition temperature Tq with increasing quark chemical potential 
is also observed on the lattice. The lattice results confirmed that the quark fluctuation in the isovector 
channel xi contrary to Xq does not exhibit a peak structure at To and shows a rather weak dependence on 
\i q as seen in Fig. 1-right. Such a properties of Xq an d Xi are expected when approaching the chiral endpoint 
with increasing fi q . However, the above behavior of Xq an d Xi can be also quantified by the regular part of the 
free energy due to the enhanced contribution of resonances in the near vicinity of the transition temperature 
|18l 1201 1251 12f>j . The above became clear in Fig. 1 where the LGT results on T and fi q dependence of 
different quark susceptibilities in confined phase are seen to be quite successfully described by the resonance 
gas partition function . The apparent agreement of the LGT results with the hadron resonance gas on the 
level of the equation of state and the susceptibilities in the temperature range T < Tq indicate that increase 
of the Xq fluctuations observed in Fig. 1 with increasing /i q at To is a necessary but not sufficient condition 
to verify the existence of the chiral endpoint. In the following we will argue in model calculations that to 
verify the appearance of the TCP in the QCD phase diagram would require a non-monotonic behavior of 
the Xq when going along the phase boundary. 

Having in mind the importance of the quark number fluctuations as a probe of phase structure of QCD 
as well as the above lattice results the main scope of this article is to consider the properties of different quark 
susceptibilities in terms of an effective chiral model. Of particular interest is the behavior of the quark number 
correlations in different channels along the boundary line and in the near vicinity of the chiral endpoint. 
The calculations will be done in terms of two-flavor Nambu-Jona-Lasinio (NJL) model formulated at 
the finite temperature and chemical potentials related with baryon number and isospin conservation. 
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Figure 2. The NJL model phase diagram in the chiral limit for Gy' = 0,0.3 and 0.6 G 3 . The dashed 
(solid) line denotes the second-order (first-order) transition line. The tricritical point indicated by a dot (•) 



is located at {T, fi q ) = (65,275) MeV for = and (T,/i,) = (42,305) MeV for Gy' = 0.3 G s . The 

results are shown for vanishing isovector chemical potential. 
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1. The net quark and isovector fluctuations in the NJL model 

The thermodynamics of the NJL model at the finite temperature and non vanishing net-quark and isospin 
chemical potentials is obtained from the partition function Z(T, [i qi fJ-i,V) formulated as a generating 
functional in the Euclidean space. In the mean field approximation jS] the partition function is obtained 
from the effective Lagrangian 



L = ${id -M + /} 7o )V> - 777- ~ ™) 2 + — 7s)(^ q ~ M 9 ) 2 + — m (fii ~ M/) 2 , (1) 

in which the thermal averages (?/>t ' 7^) have been neglected. The strength of the constituent quarks 
interactions in scalar and vector channels are parameterized by the effective coupling Gs and Gy respectively 
which have dimensions of length square. To distinguish quark-antiquark interactions in the iso-scalar and 
iso-vector channels we have defined two independent couplings Gy^ and Gy^ respectively. In Eq. Q we 
have also introduced a dynamical mass M and shifted chemical potential jl which are obtained from 

M = m- 2G S (Hj) , A = Mg + M/t 3 , (2) 

with 

fi q = fi q - 2GP (^70^) , fii = lii- 2G^ V) (^r 3 7o V) • (3) 
In the mean field approximation the thermodynamic potential of the NJL model is obtained in the following 
form 

fi(T >W M,/i) = Slf(T,ii;Mf,jif) + 

(M - m) 2 - —^(ftq - tiqf - —RKifa - M/) 2 , (4) 
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Figure 3. The quark number susceptibility \q f° r different values of net-quark chemical potentials /i q as 

a function of T/To in the chiral limit. The x? frG °' i s the quark number susceptibility for ideal quark gas 
and To = 177 MeV is the transition temperature at fi q = hj = 0. The calculations correspond to isospin 



symmetric system and the vector coupling constant G 
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'|p| 2 + Mj is a quasiparticle energy and '{T,jlf) = yl + cxp[(_E/ =f 

distribution function for the particle (+) and anti-particle (— ). 

The condensates appearing in Eqs. © _ are obtained from the conditions to minimize the 
thermodynamic potential with respect to the dynamical mass and the shifted chemical potentials, d£l/dM = 
dil/djlq = dVl/djli = 0. From the above stationary conditions and from Eqs. 101-® one calculates the 
quark condensates as the solution of the following gap equations: 



Mf =m f + 4G s iV ( 



f=u,d 



m = P-q + 
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The above gap equations together with the potential |0J arc sufficient to describe within the NJL model the 
thermodynamics and the phase structure of an effective quark medium at the finite temperature and the 
net-quark and isovector chemical potentials. 

Fig. [2 represents the phase diagram of the NJL model related with the chiral symmetry restoration 
in the (T, ii q )-pl&ne obtained in the isospin symmetric system and in the limit of vanishing current quark 
masses. The location of the boundary line that separates chirally broken from the symmetric phase was 
found from the requirement of vanishing dynamical quark mass M(T,fi q ) = when approaching from the 
side of the broken phase. Along the boundary line the order of the chiral symmetry restoration transition is 
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not unique. In the high quark density regime the phase transition is of the first order and terminates at the 
finite T and fi q as the second order transition corresponding to the position of the tricritical point (TCP). In 
the NJL model under the mean field approximation the location of the TCP is determined by the condition 
of vanishing of the second and the fourth order coefficients in the Taylor expansion of the thermodynamical 
potential: f2(M, fj, q , T) = a + a 2 M 2 + a 4 AI 4 + 0(M 6 ) applicable in the limit of M -> 0. For temperatures 
above the TCP the transition stays of the second order as expected in the Ginzburg-Landau theory. 

The position of the boundary line and the TCP clearly depends on the values of the model parameters 
[281 1291 13171 13 1| . In Fig. [3 the critical line was calculated for different strength of the vector coupling Gy 
at fixed values of Gs and the momentum cut-off A. As can be seen in Fig. [21 an increase of Gy results a 
decrease of T c and a shift of fi q towards larger values. This is to be expected [35] as the vector couplings 

Gy ' and Gy are related with repulsive interactions of constituent quarks. The position of the TCP is also 
shifted to the lower temperature and higher [i q with increasing coupling in the vector channel. A similar 
modification of position of the phase diagram and the TCP is observed with change of the scalar coupling 
Gs and the cut-off A. It is interesting to see that for a sufficiently large Gy the TCP disappears from the 
phase diagram and in the whole parameter range the phase transition is of the second order. 

1.0.1. Quark susceptibilities near the phase boundary With the thermodynamic potential and the self- 
consistent gap equations one can calculate the net-quark number and isovector susceptibilities and study 
their sensitivity and behavior near the phase transition. 

To calculate Xq an d Xi from Eq. © we have to take into account that the dynamical masses Mf and the 
shifted chemical potentials jlf are implicitly dependent on [i q , fii and T. Consequently the susceptibilities 
Xq i are controlled by derivatives of Mf and jlf 




where we have suppressed for simplicity the T and jlf dependence of rif distributions. 

The NJL model does not exhibit confinement properties of QCD. Thus, there are no hadronic bound 
states and resonances in the chirally broken phase in the NJL medium. Instead, we are dealing with 
constituent quarks which can be viewed as quasi-particles with the temperature and density dependent 
mass. In the chirally symmetric phase composition of the medium in the NJL model is not changed. At the 
chiral transition the dynamical quark masses Mf vanish and above T c the medium is populated by interacting 
massless quarks. In addition, due to the momentum cut-off there is a suppression of large momentum quark 
modes which is particularly efficient at high temperature. The differences in the mass spectrum of the 
NJL model and the QCD as well as suppression of the particle thermal phase-space will result in different 
quantitative properties of quark number fluctuations. However, this does not exclude some possible common 
features of susceptibilities in QCD and in the NJL model related with the restoration of the chiral symmetry. 

Fig. shows the quark number susceptibility \q as a function of T for different values of [i q . The 

calculations correspond to the chiral limit and were done for the fixed value of the vector coupling 

(s) 

constant G v = 0.3Gg. The net-quark fluctuations are normalized to that one for an ideal quark gas 
x ^=N c N f (Ty3 + ^). 
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Figure 4. The quark number (left) and isovector (right) susceptibilities Xq,l along the phase boundary 
line as a function of the transition temperature T c . In the left-hand figure the solid line denotes the \ q 
approaching from the chiral broken phase and the dashed line from the symmetric phase. The vertical 
dotted-line indicates the position of the tricritical point. The calculations were done in the chiral limit in 
the isospin symmetric system with the vector coupling constant Gy = 0.3 Gs. 



There are generic features in the temperature dependence of Xq f° r different values of fi q . From Fig.[3]it 
is clear that the model exhibits a phase transition for all values of [i q . The transition temperature is strongly 
/z g -dependcnt and decreases with increasing fi q . However, there is an essential difference between the critical 
behavior of Xq a t vanishing and at the finite fi q . For \x q ^ the susceptibility exhibits discontinuity at T c 
which increases with increasing fi q . For [i q = such a discontinuous structure disappears and instead at T c 
the susceptibilities shows a non-analytic structure which results as discontinuity in higher moments of net- 
quark number fluctuations. Such a properties of Xq are consistent with that expected if the phase transition 
is of second order and belongs to universality class of three-dimensional 0(4) symmetric spin models |21II25| . 
To see it, let us construct an effective thermodynamic potential in the near vicinity to the chiral transition. 
The relevant field is here a constituent quark that carries the dynamical mass M. Performing the Taylor 
expansion of f2(M, T, /j, q ) as power series around M ~ one gets: 

Sl(T,tt q ,M) ~ n (T, N ) + ^a(T,n q )M 2 + h(T, N )M 4 . (11) 

Thus, f2 has a structure of the Ginzburg-Landau (GL) potential where the effective quark mass acts as the 
sigma field. Following the GL theory and applying the mean field approximation the effective potential (If f ll 
is rewritten as 

n(T^ q ,M ) = n (T, N ) \ a l£^ ■ (") 

where we have used the stationary condition dfl/dM\M = and introduced the stationary point 

Mo = y/—a/b. Approaching the critical line from the symmetric phase with M = 0, the quark number 

susceptibility Xq 1S obtained from Eq. lfTT|) as 

(sym) = d 2 ^o . . 

In the GL theory the second-order phase transition line (the 0(4) critical line) is determined by the 
requirement that the coefficient a = and b ^ in Eq. (|llfl . In the near vicinity to the critical point (T c , \i c q ) 
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the coefficient a(T, can be parameterized in the MF approximation as 

a(T lN )~C(T-T c ) + D( N -^ q ), (14) 

where C and D are independent of T and /i q . Approaching the 0(4) critical line from the broken phase, the 
quark susceptibility is calculated from Eqs. I|12fl and (|14fl as: 

(broken) = d 2 »o _ D 2 , ■. 

X « 8,4 2b(T, N )- { > 

Comparing Eqs. (|13[l and l|15|) it is clear that the second term in (|15|l gives just a discontinuity of Xq across 
the 0(4) critical line at the finite fi q . While at fx q = the coefficient D = and x (sym) = x (brokon) at T = T c . 

Considering the phase diagram in Fig. [2] we have already indicated that for large \i q the NJL model 
experiences the TCP. The right pannel of Fig. shows the properties of the quark number susceptibility in 
the near vicinity and at TCP. In the GL theory the position of the TCP is characterized by the condition 
of vanishing a(T,/i q ) and b(T,fi q ) coefficients in the effective potential Q11JI. Consequently, from Eq. i|15|) it 
is clear that the quark fluctuations Xq should diverge at TCP. These expected properties of Xq are clearly 
seen in Fig. [21 With the present choice of parameters the TCP is located at (T c , [i q ) — (42, 305) MeV where 
X q — * oo. Crossing TCP towards larger /j, q the second order 0(4) chiral phase transition is converted to the 
first order where Xq is finite and has a gap at the critical temperature. 

From the perspective of heavy ion experiments the properties of different susceptibilities are of particular 
interest. This is because, fluctuations related with the conserved charges are experimentally directly 
accessible . Since these are also observables that are sensitive to the critical behavior, knowledge of the 
properties of susceptibilities along the phase boundary line could give insights how to verify the QCD phase 
transition experimentally. Clearly, the quantitative structure of the phase diagram and the position of the 
chiral endpoint is model dependent. Thus, also the position of the QCD boundary line could be very different 
than that found in the NJL model. However, the model study could still answer a phcnomcnologically relevant 
question how to observe and how large is the critical region along the phase transition where the fluctuations 
are sensitive to the singular structure at the tricritical or chiral endpoint. 

Fig. 01 shows the net-quark susceptibility Xq along the phase boundary line from Fig. |2 The Xq is 
quantified as a function of the chiral phase transition temperature T c . The appearance of the TCP in the 
phase diagram results in a non-monotonic behavior of Xq along the transition line. There is a window 
of AT C ~ 30 MeV above and Afj, q ~ 10 MeV around TCP where the fluctuations are sensitive to the 
appearance of the tricritical point. If the TCP was absent in the phase diagram then the Xq would be a 
monotonic function of T c along the phase boundary as shown by a dashed-dotted line in Fig. Such a 
behavior is seen in Fig. 4-right in the isovector susceptibility which is not sensitive to the appearance of 
the TCP in the phase diagram j33j . The non-monotonic behavior of Xq is onr y seen from the side of the 
chirally broken phase. Approaching T c from the chirally symmetric phase results in continuous behavior 
of Xq along the boundary line. This is because the xi^™^ is finite at the chiral transition in the whole 
parameter range as seen in Fig. A difference between Xq Sym ^ and ^ brokon ) calculated along the phase 
transition line measures the magnitude of discontinuity at the phase transition. Fig. 0] shows that at [f q = 
this discontinuity vanishes and at T c = would be the largest if the phase diagram did not experience the 
TCP. 

In nucleus-nucleus collisions any change of the temperature and chemical potential is correlated with 
the corresponding change of the c.m.s collision energy yfs |34j . An increase of y/s results in increasing the 
temperature T and decreasing the quark chemical potential fj, q . Thus, the critical region around tricritical or 
chiral endpoint (A/Lt°, AT C ) can be converted to that in the c.m.s. -y/s-energy in A-A collisions. Admitting 
that the relation of fi q and T c with y/s are the same as for chemical freezeout parameters extracted from 
experimental data |34) the AT C ~ 30 MeV window around the TCP would correspond to As ~ 1 A-GcV. 
Thus, within our crude estimates, we can expect that to observe any remnant of critical fluctuations in A-A 
collisions one would need the -y/s energy step to be within a range of 1 AGcV. 



Charge Fluctuations along the QCD phase boundary 



8 



We have discussed the properties of quark fluctuations in terms of Nambu-Jona-Lasinio (NJL) model. 
The model was formulated at finite temperature and chemical potentials related with the conservation 
of baryon number and isospin. Applying the mean field approach, we have shown how the fluctuations 
of different quark flavors are changing across the phase boundary. Such a study is interesting from the 
perspective of heavy ion phenomenology and the lattice gauge theory. In the first case we have indicated and 
quantified the non-monotonic structure of the net-quark number susceptibility along the phase boundary as 
the method to identify the position of the chiral endpoint. We have also discussed the critical region around 
this point in the context of heavy ion phenomenology. 
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discussions with H. Fujii and B. J. Schacfer. K.R. acknowledges stimulating discussions with S. Ejiri, 
K. Rajagopal, E. V. Shuryak and M. A. Stephanov. The work of B.F. and C.S. were supported in part 
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